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D ' Abstract 

We consider flat Friedmann-Lemaitre-Robertson- Walker cosmological models in the 
framework of general scalar-tensor theories of gravity with arbitrary coupling functions, 
set in the Jordan frame, in the cosmological epoch when the energy density of the ordinary 
dust matter dominates over the energy density of the scalar potential. Motivated by 
Qi cosmological observations, we apply an approximation scheme in the regime close to the 

' so-called limit of general relativity. The ensuing nonlinear approximate equations for the 

scalar field and the Hubble parameter can be solved analytically in cosmological time. 
This allows us to distinguish the theories with solutions that asymptotically converge to 
general relativity and draw some implications about the cosmological dynamics near this 

00 ' limit - 

O 
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1 Introduction 

Various cosmological observations of our Universe can be fairly well accommodated within the 
ACDM concordance model pQ based on the theory of general relativity (GR). However, there 
is still a number of viable alternative theories which also manage to conform sufficiently well 
with observational data [2]. One such family of theories is provided by scalar-tensor gravity 
(STG) [3] where gravitational interaction is mediated by an extra scalar degree of freedom 
\1/ in addition to the usual tensor ones. In the so-called Jordan frame and Brans-Dicke like 
parametrization an STG is characterized by two arbitrary functions, the coupling function 
w(\I>) and the scalar potential V(ty). As has been discussed by many authors previously pfllH], 
for a range of choices of u> and V the cosmological evolution of dust and potential dominated 
STG models naturally converges close to the one expected from GR. Yet, at the same time 
STG models may also offer a possibility to explain small observational differences from pure 
GR ACDM behavior, e.g. the possibly variable effective barotropic index of dark energy [6] as 
hinted by some observational data [7], deviations in the growth of perturbations [8J, etc. 
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The aim of the current paper is to narrow down the class of STG models that can lead to 
observationally viable cosmologies (i.e., spontaneously evolve close to GR), and by explicitly 
finding the general solutions applicable in this regime to provide a basis for further direct 
checks with observational data. It is a follow-up work to our recent papers pi [lOj [H] where 
we investigated Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmological models in the 
framework of general STG with arbitrary coupling function and scalar potential in the era 
when the energy density of the scalar potential dominates over the energy density of ordinary 
matter. There we presented, justified and applied an approximation scheme for the scalar field 
equation to capture the scalar field dynamics near the GR limit. In the present paper we 
supplement these studies with analogous investigations for the the cosmological epoch when 
the energy density of the ordinary dust matter dominates over the energy density of the scalar 
potential. The presence of an extra dynamical quantity (matter) in the system makes the 
procedure now a bit more complicated, yielding two non-linear equations which explicitly 
contain time but which can be nevertheless solved analytically. In a realistic cosmological 
scenario the dust dominated epoch should be patched together with the potential dominated 
era (as well as with an account of the early universe). For some related recent studies see Refs. 

In section 2 we recall STG FLRW equations. In section 3 we motivate and apply to the 
dust matter dominated era the approximation method worked out in Ref. [9]. The resulting 
nonlinear equations are solved analytically in cosmological time in Sect. 4. Comparison with 
earlier results and implications for selecting a model of STG viable in cosmology are discussed 
in section 5. Finally, section 6 provides a summary and a brief outlook. 



2 The equations of scalar-tensor cosmology 

We consider a general scalar-tensor theory in the Jordan frame given by the action functional 

+ S m {g ta ,,Xm)- (1) 

Here ufif) is a coupling function, V M denotes the covariant derivative with respect to the 
metric g^ u , k, 2 is the non- variable part of the gravitational constant, and S m is the matter 
contribution to the action as all other fields are included in \ m . In order to keep the effective 

2 

gravitational constant 871G = 4r positive, we assume that < ^ < 00. 

The field equations for the Friedmann-Lemaitre-Robertson- Walker (FLRW) line element 

ds 2 = -dt 2 + a(t) 2 ( + r 2 (d6 2 + sin 2 8d<p 2 )) (2) 
\ 1 — kr z J 

with curvature parameter k = (flat) and perfect barotropic fluid matter, p = wp, w = const., 



VR{g) - ^^V^Vp^ - 2k 2 V(V) 
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read 



\j/ 1 \j/ 2 \Jf K 2 K 2 

2H + 3H* = '2H----um----, P + -Vm, (4) 
$ = -3ff^ + ~ A(g)(2t«;(*) + 3) fr 2 + £ — (1 - 3w)p 



2k 



2 



2w(^) + 3 



2V(tf) - * 



(5) 



where H = a/a and we have introduced the notation 



for later convenience. The matter conservation law is the usual 

p + 3H (w + 1) p = 0, (7) 

it is reasonable to assume positive matter density, p > 0. 

The Hubble parameter H can be expressed as a function of \& by solving the Friedmann 
equation (JHJ) algebraically, 



g = --i±. U)^/^^)), ( 8) 

2# V 12^ 2 3^ v ; 

For later argument notice that in the limit i 2 d(1)+3) — ^ ^> ^ 7^ the system faces a spacetime 
curvature singularity, since if diverges. Only as long as (2cu(\E') + 3) l I /2 is finite are the solutions 
singularity free. 

Let us take the regime where the dominating contribution to cosmological energy density 
is provided by dust matter (w = 0) and the scalar potential can be neglected in the equations. 
The system (J3])-([TJ) is characterized by three variables H, p}, but one of them is algebraically 
related to the others via the Friedmann equation ([3]). Eliminating p yields two equations 

« = _ 3H y + l -{2u + Z)A{^ 2 + j^—j ^H 2 + 3H* - |L | , (9) 

3 \i> \i> 2 i \i/ 2 

H = --H 2 + H -u- -(2u + 3)Am — 

2 2^ 4^ 2 4 V ; v ; V 

1 / vjr ^ 2 

ZH 2 + ZH--—u> , (10) 



2(2w + 3) I * 2^ 2 

which provide the basis for the present study. 
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3 Approximate equations 



In fact, not all possible solutions of Eqs. (19|)-( TT0|) are of immediate physical interest, since 
cosmological observations give a clear preference towards a certain corner in the solutions 
space. Analysis of the anisotropies of the cosmic microwave background (CMB) radiation sets 
a limit on the variation of the gravitational constant from the recombination process till now, 
IGrcc-Gnowl < 5x 1CT 2 [131 HU , which for scalar-tensor gravity translates into 4? -C 1. In addition, 
the best fit of the CMB data indicates that at the time of recombination „ ^t, 1 ^ , q < 7 x 10 -2 

^V W rcc )+A 

[14"] . while the value today from the PPN data is bounded as 2oj ^ 1 ^ +3 < 7 x 10~ 4 [13] . 

Therefore it makes sense while considering the dust dominated cosmological era to focus 
upon the solutions near the limit (a) ^ 2 ui(l)+3) ~^ ® an< ^ O 3 ) — y ®- This assumption is 
consistent with the equations, as one can check in Eq. (J3J how the conditions (a) and (b) keep 
\1/ negligible, thus allowing \l/ to remain negligible as well. One may also argue that if changes 
in \I/ are sufficiently small, does not change dramatically and the regime expected from 

the solutions is sufficiently stable to merit investigation. 

So, let us define ^ by 

and focus upon the solutions near this point, 

V(t) = x(t) , H(t) = H*(t) + h(t) , (12) 



where H*(t) is the Hubble parameter corresponding to the cosmological evolution with \l/ 
while x(t) and h(t) are small deviations. It follows from (JT2"]) that = x(t), where we 

expect x(t) to be also small due to (b). Under the two additional mathematical assumptions, 
(c) A+ = A(ty+) and (d) ^—3 is differentiable at we can expand in series 

1 1 

+ A*x + ... mA+x, (13) 



2u;(#) + 3 2w(**) + 3 

(2oo(*) + 3)4> 2 = — -f— = *L(l + (x)) n j?- . (14) 
+ A*x + • • • A+x A*x 

The latter result actually informs us that in order to avoid a spacetime singularity — must 
not diverge, hence we should treat x{t) and x(t) as the same order (small) quantities, cf. the 
remark after Eq. (jSJ) . In passing let us remark that in our previous papers [9j [TUJ E] we have 
tentatively called (a)-(d) 'the limit of general relativity' since under these conditions the set of 
STG cosmological equations dBJ)-© reduces to those of pure GR (with a cosmological constant 
ifV(tf*)^0). 

Subjecting the H equation (ITU]) to the approximation ( TT2l gives in the first order 

H* + h = -\hI - ZHJt - ( 1 + tt-tVI — + 7n~ H ** ~ l^Hlx . (15) 
2 4** V 2A^/ x 2^ 2 

Taking the limit where the deviations x, x, h, h vanish, we define 

= (16) 
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which is familiar from the Friedmann solution of the dust dominated pure GR. It determines 
the time evolution of H* to be 

H ' = W^7Y (17) 

Here t s is a constant of integration which fixes the beginning of time scale; in what follows we 
choose t s = 0, t > 0. For late times when is finite, Eq. ffT5l) now assures that h is also 
small, at least on a par with h. To sum up, the appoximate first order equations read 



x 



2 



x = — - 2>H±x + ZAJU+Hlx , ill 
If 1 \ x 2 1 



h + 3HJl ~w, \ 1+ 2ajtj ~ + m H - x - 2 AMx (19) 

with H* given by Eq. (JT71) . Notice that due to the Eqs. (TTBj) . ( JT9l) depend explicitly on 
time t. This means that the corresponding system of first order equations is not autonomous 
and the standard phase space analysis is not applicable. However, we can straightforwardly 
integrate Eqs. (j!8p . f|T9|) in cosmological time and analyse the behaviour of solutions in the 
neighbourhood of the limit of general relativity. 

For later reference let us note that the expansion ( jT2l) can be also applied for the effective 
barotropic index, 

w eff = -1 - « ~ (h + 3H*h) . (20) 



Thus once h(t) is found, it can be plugged into the equation above to reveal how w e ff evolves 

G 
G 



in cosmological time. In an analogous manner one may also deal with § and other relevant 



quantities. 



4 Solutions in the cosmological time 

Despite its nonlinear and nonautonomous structure, one can solve Eq. f|T8|) analytically. It 
turns out that the type of the solution x(t) depends on the constant 

D = 1 + j^tf* (21) 

which characterizes the underlying STG. Then one can plug in x(t) into Eq. ffT§j) and solve the 
latter for h(t), which also yields an analytic result. Having found x(t) and h(t), it is possible 
to determine the evolution of the effective barotropic index w e ff from Eq. (120]) . and other 
quantities of interest. 

4.1 Polynomial solutions 

In the case D > the solution of Eq. (fTBj) reads 

± x(t) = - - M 2 t-^y . (22) 
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Here and below the "±" follows from an obvious invariance property of Eq. (1181) under reflec- 
tion x —x, i.e. there are solutions which lie in the regions ^ ^ (x ^ 0), respectively. 
The constants of integration M 1; M 2 are related to the initial data x* = x(t*), x* = x{t*) at 
some arbitrary time t* as 



Mi 
M 2 



+ x*(l + \/D) i(i_ 
2x/D x /±xZ 

+ x*(l - VD) |(i +v ^d) 



Now we can integrate Eq. (fT9l) to obtain 

± fc(t) = — Ml (-ay/D + &) + M| (av 7 /) + &) r^ 5 + K 



(23) 
(24) 

(25) 



Here K is another constant of integration and we have introduced constants a, b which char- 
acterize the underlying STG, 



3 + 6A*% 



3 + lOA^t 



As a result, the full Hubble parameter in the approximation under consideration is 
H{t) = j t jl ± ~ [M 2 (-av 7 ^ + &) ^ + M 2 2 (a^D + 6) + K 



The effective barotropic index reads 



± Weff(t) 



M 2 (-av^D + b\ t^ 5 - M\ (ay/D + b\ t~ 



(26) 



(27) 



(28) 



We can get a better feel of these solutions by considering their behavior at certain limits 
and points. Asymptotically at t — > oo the solutions exhibit two distinct behaviors. For STGs 
with \/D < 1 (i.e. A*\V* < 0) all cosmological solutions irrespective of their initial conditions 
monotonically approach the general relativistic dust matter FLRW cosmology, ^(i) — > = 
const., H(t) —> H*(t) = 2/(3i), w e fr(t) — >■ 0, since all first order corrections vanish at this 
limit. On the other hand STGs with \/D > 1 (i.e. A+\V± > 0) allow only solutions that will 
diverge, x{t) — > oo, h(t) — > oo, w e s(t) — > oo, meaning that solutions in these theories can 
linger near general relativity only for a certain period, while as time evolves they will leave 
and the approximation scheme will break down eventually. (The case \TD = 1 would imply 
A+ = or \V± = 0, which contradicts the assumptions (c) or < \V < oo of the present study.) 

Taking t — > the quantities x(t) and h(t) diverge for all integration constants and param- 
eters of the theory except for the special M 2 = 0, \/D > 1 case. This indicates that generally 
the solutions can not start near the limit of general relativity and only dynamical evolution 
can bring them close to it. 

The solution (T22]) also informs us that if the integration constants Mi and M 2 are both 
positive or both negative, then at a finite moment 



M 

Mi 



1 



T±\ > 



(29) 
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the corresponding solutions can go through 



± x(t b ) = , ±x(t b ) = , ±x(t b ) = 2DM 1 M 2 t b - i > . (30) 

At t b these solutions do not stop at x — 0, but bounce back, i.e. the solutions coming from 
the x < region return to the x < and similarly the solutions coming from the x > region 
return to the x > region. There is no crossing from x < to x > or vice versa. In terms 
of the initial data at some arbitrary time the bouncing solutions satisfy 

^iA<±i t (l-\/^), or ± x*t* < =Fx*(l + Vd) , (31) 

as can be inferred from Eqs. 
In addition, at 



i 



_ M 2 (1 + ^D) 



the solutions may pass through 

± x(t c ) = , ±x(t c ) = , ±x(t c ) = -2DM 1 M 2 tf . (33) 

This happens for two types of solutions: if sig^A^) = sign(M 2 ) for a/Z) < 1 and if sign(Mi) ^ 
sign(M 2 ) for \J~D > 1. The first type encompasses all \/D < 1 solutions which at t b have 
bounced back from x = 0, now at t c > t b they turn around again to proceed asymptotically 
towards x = 0. The second type comprises of the \/D > 1 solutions which never get to x — 0, 
the moment t c marks their closest reach to x = before starting to flow away. 

Therefore, in summary, the following picture emerges. If VD < 1 the solutions with initial 
conditions ( 13T1) first apporach x = 0, then at t b reach x = and bounce back, further at t c 
turn towards x = again, to get there asymptotically as t — > oo. The a/D < 1 solutions with 
initial conditions outside the ranges given by ( 13~T|) converge to x = monotonically. If a/D > 1 
the solutions with initial conditions ( l3Tj) initially approach x = 0, then at tf, reach x = to 
bounce back and flow away. The a/D > 1 solutions with initial conditions outside the ranges 
given by (13TT) move towards x = 0, but before reaching it turn around at t c and leave. An 
exceptional case is the \[D > 1, M 2 = solution which starts at x = and monotonically 
flows away from it. 

As is evident from f l25|) and fl28l) the behavior of h{t) and w e fj (£) is not synchronous with x(t). 
However, one can make some simple generic observations taking into account the definitions 
( 12T1) . ( )26|) and the basic assumption < ^ < oo. First, as t — > the quantity w e fr(t) — > — oo 
if | < D < 1, while w c fj(t) — > +oo if0<D< = orZ)>l. Also, as we have noted above, 
t — > oo takes w e ff (t) — )■ if D < 1, while w c g(t) — >■ +oo if D > 1. Further, there can be specific 
moments 

/ Mf (a/Do + 6) \ 575 . . , , , 

«- = ^ ( _^ fl + 6) j ' Wcff = ° ' Wcff(td) + ° (M) 

M| (1 + VD)(VDa + b) 



and 



(1 - y/D)(-y/Da + b) 



w eff (t e )^0, w cff (t e ) = 0. (35) 
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One can check that t d > if < D < |, while t e > if < D < | or D > 1. 

Thus, the picture is the following. For < D < | the barotropic index w e ff(t) approaches 
w e g = from above, and passing this value at t d , then later at t e the quantity w e s(t) starts 
to increase again, and will asymptotically converge to the vanishing value. For | < D < 1 
the solutions exhibit a monotonic growth for w c g which closes in to the w e g = dust matter 
regime from below. The generic D > 1 solutions start with decreasing w c g, which reaches 
its lowest (and positive) value at t e , but after that w e g- starts to increase again. The D > 1, 
M 2 = solution is an exception, here w e g starts from and keeps increasing in time. 

4.2 Logarithmic solutions 

In the case D = (A+V* = -§) the solutions of Eqs. (ITSI-CT read 

±x(t) = - (Mi hit- M 2 ) 2 , (36) 



t 



±h(t) 



Mi 



^-(\nt) 2 + (M! - M 2 )lnt + K 



(37) 



±w cff (t) = -^(Milnt + Mi-M 2 ) • (38) 

Here Mi,M 2 as well as K are constants of integration, fixed by the initial data = £*(£*), 
(U) at some arbitrary time t* as 

Mi = ^5^*1, (39) 
~ + a;*) In — 2x* ~| 

M2 - (40) 

Unless Mi = the solutions exhibit the same generic behavior. They start by approaching 
x = 0, at 

M 2 

t b = e^h (41) 

reach 

± x (t b ) = , ±x(t b ) = , ±x(t b ) = 2M^t b 3 > (42) 
and bounce back, later at 

2A/ 1 +A"/ 2 

t c = e ~^~ > t b (43) 

the solutions pass through 

±x{t c ) = AMli-\ ±x(t c ) = 0, ±x(t c ) = —2M\t~ z , (44) 

and return flowing towards x = reaching it asymptotically in time. The effective barotropic 
index starts by decreasing from a positive value, experiences 

j£l-jf2 _ _ M? ~ r, 

t d = e , w cff (t d ) = 0, w cS (t d ) = (45) 
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but then starts to incease again at 

t e = e"i = t b , w eff (t e ) = -i^-i- 1 , w cff (t e ) = , w eff (f e ) = ^-tf (46) 

reaching w e g- = from below asymptotically in time. In the Mi = case the solutions approach 
x = monotonically while w e g is always zero. 



4.3 Oscillating solutions 

In the case D < solutions of Eq. ( JT8l) read 

±x(t) = 



.\\ sin ( -y^D|lntj -iV 2 cos QA/]D|lnt 



n 2 



(47) 



where the constants of integration Ni and N 2 determined by the initial conditions X* = x(t*), 
(U) at some arbitrary time as 



N 2 



\D\^/±xZ 



Vl^lv 7 ^ 
The corresponding solution of Eq. (|T9|) is 



|x*t* + £*) cos ( -A/FDjlnt*J + x*a/|-D| sin ^-a/|D| hit* 
4A + x.) sin ( i v^fln*,) - v^cos (\ yffilat. 



(48) 
(49) 



±h(t) 



2 

3t2 



(iV 2 2 - JV?) 



£u/|Z>| 



- AqiV 2 fej sin (v^Djlnt 



+ (V 2 2 - Nl)^ + A^iV, av^Dj^) cos (y/\D\ lnt) + 



(50) 



where is another constant of integration and a, 6 are given by Eq. (126]) . The full Hubble 
parameter now consists of small oscillations around the GR FLRW dust cosmological model 



Hit) = -\l±- 
W 3t I t 



{Nl - Nl) 



j^i - N X N 2 bj sin (y\D~\lnt 
+ (V 2 2 - iV 2 )^ + av^Dj^) cos (y/\D\ 



lnt) +K 



(51) 



The effective barotropic index reads 



± Weff(t) 



^(Af - A\ 2 )^ + N^a^/lJT^ sin (\/^D|ln 



+ ( (Af - Nl f^ Dl - NiN 2 b ) cos (y\D\\nt 



(52) 



The behavior of all these solutions is fairly simple as they approach the general relativistic 
dust matter cosmology in the manner of damped oscillations. At the moments 



t b = exp 



\D\ 



arctan 



N2 



+ 



2nn 
^f\D\ 



(53) 



the deviation x(t) of the scalar field passes through 

±x(t b ) = 0, ±x(t b ) = 0, ±x(t b ) 
i.e., bounces back from x — 0, while at the moments 



—D 



{Nl + iV 2 2 )t- 3 > 



t c = exp 



\D\ 



: arctan 



N 2 + N iy f\D\ 
N 1 -N 2 J\D\ 



2nir 

\Ad\ 



(54) 



(55) 



it passes through 
± x(t c ) 



1 — D 



{Nl + N 2 2 )r c \ ±x(t c ) = 



±x(t c 



D 



(N* + N%)t: 3 <0, (56) 



i.e., turns around and evolves towards x = again. The amplitude of the deviations monoton- 
ically decreases while the period monotonically increases. The behavior of w e s(t) is analogous, 
but not synchronous with x(t). It is characterized by 



Weff(td) 



exp 



\D\ 



arctan 



(iVf - N?)ay/\D\ - 2N x N 2 b 
(Nl - Nfjb + 2N l N 2 a^/\D\ 



nir 



\D\ 



0, w eff (t d )^0 



(57) 



and 



exp 



\D\ 



arctan 



(Nl - Nf)(b + a)D + 2N 1 N 2 (b + aD)y/\D\ 
(Nl - Nl)(b + Da)^\D\- 2N x N 2 (b + a)D 



+ 



nit 



\D\ 



±W e ff(t e ) 

where 



-D(Nl + Nt)VW^D , 
S _ = * e , w e fr(t e ) - U 



q = (-l) n sign ((Nl - N*)(b + Da)y/\D\- 2N Y N 2 (b + a)D 



(58) 



(59) 



meaning oscillations around w e ff = with exponentially decreasing amplitudes and expone- 
tially increasing period. 



5 Discussion 

5.1 Comparison with earlier results 

STG dust cosmology equations near the GR limit were investigated several years ago in the 
Einstein frame by Damour and Nordtvedt [4J. By invoking an analogy with a mechanical 
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particle with time- dependent mass, they demonstrated that in the case of coupling function 
(2a; (\&) + 3) -1 / 2 = a(cp) = kip, k = const, the type of a solution for the Einstein frame scalar 
field tp(p) with the evolution parameter p = (2/3) In t depends on the numerical value of the 
model-dependent constant k: the solution is exponential in time parameter p, i.e. polynomial 
in cosmological time t if < k < 3/8, linear-exponential if k = 3/8 and oscillating if k > 3/8. 

In our earlier papers [9] we investigated the Jordan frame scalar field equation close to the 
GR limit in the linearized approximation, found the fixed points and calculated the eigenvalues 
which determine the type of solutions around these fixed points. Our results were qualitatively 
similar to those of Damour and Nordtvedt [4], but the critical value of the model-dependent 
parameter turned out to be 3/16 instead of 3/8. 

In the present paper we refined the analysis and found solutions in the nonlinear approx- 
imation for the Jordan frame scalar field \l/(t) in the cosmological time t and obtained the 
critical value of the model-dependent parameter to be given by A±^+ = —3/8. It is in ex- 
act agreement with the results of Damour and Nordtvedt, as the transformation between the 
Einstein and the Jordan frame quantities 

{dvr= ^m±i m * (60) 



gives 



k ^ a I 
dip ' 



2^ du 



-4**. (61) 



_(2w + 3) 2 c^_ 

It follows that the approximation used by Damour and Nordtvedt |4j in the Einstein frame is 
congruent with our nonlinear approximation in the Jordan frame and thus can be considered 
as an additional justification for our expansions ( 1T2|) . ( ITBl . 



5.2 Combining the dust and potential dominated eras 

In the present paper we focussed upon the dust dominated cosmological epoch in the frame- 
work of STG with negligible scalar potential. In principle, this epoch could be followed by 
a scalar potential dominated epoch with insignificant matter density that we investigated by 
similar methods in our earlier papers jTDJ [11] . In both cases we assumed that the cosmological 
model has evolved towards the GR point ^ ( II II) since this is strongly indicated by different 
contemporary observations, and we solved field equations in a nonlinear approximation in the 
neighbourhood of this point. Let us now combine the conditions on the parameters of the 
models with the aim to view different epochs as parts of a single cosmological scenario. 

In both cases there are general conditions for solutions to converge towards the GR value 
asymptotically in time: in the dust dominated model it reads (see Sec. 4) 



d 



d^ \2o;(#) + 3 
and in the potential dominated model [10] 



< 



(62) 



>0, 



* dV 
2VdV 



< 1 



(63) 
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Epoch 



Solutions 



p dominates 



oscillating 



logarithmic 



polynomial 



D < 



£1 = 



< D < 1 



V dominates 



oscillating 



linear — exponential 



exponential 



B < 



B = 



< B 



Table 1: Classification of the qualitative behavior of solutions of the scalar field and cosmologi- 
cal expansion while converging to the GR limit in the dust matter (p) dominated and potential 
(V) dominated epochs, determined by the parameters D (I2TI) and B ( )64|) . which characterize 
the underlying STG. 

The converging solutions can be classified according to the numerical value of a model- 
dependent parameter as summarized in Table 1. As discussed in Sec. 4, in the dust dominated 
epoch the behavior of the scalar field is determined by the quantity D (12"T|) characterizing 
the STG model: the solutions are oscillating if D < (A±^+ < —3/8), logarithmic if D=0 
(A*ty* = -3/8) and polynomial if < D < 1 (-3/8 < A+ty* < 0). In the scalar potential 
dominated models the corresponding classification can be given in terms of a model-dependent 
quantity 



as follows: the solutions are oscillating in cosmological time if B < 0, linear-exponential if 
B = and exponential if B > [llj. The same behavior carries over to the cosmological 
expansion as encoded in the Hubble parameter H or barotropic index w e g, i.e. polynomial, 
oscillating etc. convergence towards the dust FLRW values in the matter dominated epoch or 
de Sitter values in the potential dominated epoch, correspondingly. 

A realistic STG cosmological scenario compatible with observations would better need to 
have GR as an attractor in both dust dominated and matter dominated regimes. Therefore 
for a credible STG both conditions ( 1621) and ( 1631) must be satisfied, thus constraining the set 
of functions oufi?) and V(\&) one can consider for constructing a viable model. The next filter 
is provided by qualitatively different behaviors among this converging class of models, e.g. 
depending on oj(ty) and V(\&) the evolution may be oscillating in the dust dominated and 
exponential in the potential dominated epoch, etc., which might be possible to detect in future 
observations. 

6 Summary and outlook 

In this paper we have considered generic Jordan frame STG flat FLRW cosmological models 
in the dust dominated era with negligible scalar potential near the limit of general relativity 
as favored by various observational constraints. We derived and solved nonlinear approximate 
equations for small deviations of the scalar field and cosmological expansion from their GR 
limit values. Depending on the scalar field coupling function u(*f?) the models fall into two 




(64) 
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classes where either all solutions approach GR asymptotically in time or only a single fine-tuned 
solution does. The models with universally converging solutions come in three characteristic 
types: polynomial convergence, logarithmic convergence, and damped oscillations around gen- 
eral relativity. 

The approximation scheme assumes that the first derivative of 2u j(^)+3 w - r -t- ^ evaluated at 
the GR limit (II lj) is nonvanishing and finite, while the higher dervatives do not diverge. Then 
the only parameter characterizing the underlying distinct STG which enters the approximation 
equations and the analytic solutions is the value of the first derivative. Thus in principle the 
present study encompasses a very large generic family of STG models and in this sense has 
wider applicability than considering example models, equvalent to a particular form of u(^/) 
chosen. 

The class of STGs where the GR limit is an attractor for the nearby solutions is of interest 
because there is a dynamical mechanism naturally driving the solutions to satisfy observational 
constraints. So, combining the results of the present work on the dust dominated epoch with 
earlier results on the potential dominated regime, provides a reasonable viability filter for STG 
models in terms of the conditions (IB"2"j) . ( 16"3j) . 

On the other hand the converging solutions still have their characteristic small deviations 
from the ruling ACDM scenario. Given the generic analytic solution for the cosmological 
expansion and the corresponding effective barotropic index near the GR limit, it remains as 
a future work to face it with actual data and to draw observational constraints on the STG 
models. Similarly, the expansion history enters as background evolution in the equations for 
the growth of perturbations, which leads to another line of investigation. Finally, the dust and 
late-time potential dominated epochs must be patched together with an account of the early 
universe. 
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